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. The mid-infrared spin-wave spectrum of antiferromagnetic YBa2Cu306.o 

Q>^ I was determined by infrared transmission and reflection measurements {k\\c) 

' at r=10K. Excitation of single magnons of the optical branch was observed 

Ch ■ at Eop = 178.0 meV. Two further peaks at 346 meV (ss 1.94 E^p) and 470 meV 

■ (~ 2.6 Eop) both belong to the two-magnon spectrum. Linear spin wave theory 

. is in good agreement with the measured two-magnon spectrum, and allows 

[ to determine the exchange constant J to be about 120 meV, whereas the 

^ I intrabilayer coupling J12 is approximately 0.55 J. 

^ '. 74.72.-h, 74.25. Gz, 74.25.Ha, 78.30.-j, 75.30.Ds 

^ : 

2 I High temperature superconductors are basically layered copper-oxide materials. It is 

O I widely accepted that the relevant electronic degrees of freedom are confined to copper-oxide 
planes. The number of CUO2 planes per unit cell varies: e.g., La2-xSr3;Cu04 exists in a 
single plane form with a large spacing between planes of f»13.2A, and YBa2Cu306+x has a 
double layer structure with intra- and interbilayer spacings of f«3.3A and 8.5A, respectively. 
Electronic correlations, and hence spin dynamics , may depend on the type of stacking of 
^ [ the planes. More specifically, a sizable coupling J12 between spins on adjacent planes of a 
Q ' bilayer will influence the spin excitation spectrum as well as the nature of the ground state. 
O ■ This may have been seen already in doped compounds: the normal state spin susceptibility 
>• . of La2_a;Srj.Cu04 extrapolates to a finite value at zero temperature, whereas it extrapolates 
^ ! to zero for YBa2Cu306.6 0- This may be interpreted as a signature for the opening of 
^ [ a spin excitation gap in YBa2Cu306.6 at low temperatures |^ — a behavior certainly not 
- - ' encountered in Fermi liquids. Further, a spin density wave ordering for La2_xSr2:Cu04 has 
been proposed, but for YBa2Cu306.6 a singlet pairing of spins in adjacent CUO2 planes 
with strong antiferromagnetic fluctuations within a plane Such a scenario seems to 

require an unrealistically large J12 > 2.5 J 0, where J is the in-plane exchange coupling 
of the Heisenberg Hamiltonian supposed to describe the low energy spin dynamics of a 
single bilayer for zero doping (x = 0). However it was argued that, for finite doping, the 
itinerant carriers destroy the antiferromagnetism of the insulating phase and, therefore, 
much smaller values of J12 will produce a singlet interplane pairing in the conducting phase 
of YBa2Cu306.6- 

Up to now, no experimental evidence has been given of a sizable bilayer coupling ( J12 ~ J) . 
In neutron-scattering experiments on YBa2Cu306+x, the in-plane coupling was determined 
from the dispersion of acoustic spin- waves and was found to be extremely large (J = 120 ± 
20 meV 0, J = 150 meV 0, both for x = 0.15). Yet, no optical modes have been found 
for energies up to 60 meV WIM, suggesting a bilayer coupling of J12 > 8 meV. In Raman- 
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scattering experiments on YBa2Cu306+a; a two-magnon peak was observed pO|JTT[] . J was 
found to be consistent with the neutron-scattering data, whereas J12 was neglected. 

In this Letter we report the first observation of an optical magnon peak and of the two- 
magnon spectrum in infrared spectroscopy of YBa2Cu306.o- They allow to determine both 
J12 and J . Two resonances of the two-magnon spectrum confirm the location of the single 
magnon peak. Further, understanding the 'antiferromagnetic limit' YBa2Cu306.o will be 
crucial to interpret the excitations of carriers in doped YBa2Cu306+x [0]- 

Due to the high resolution and the wide spectral range, optical spectroscopy is a powerful 
method to determine precisely energies of spin- waves. But, compared to other optical exci- 
tations such as infrared-active phonons, intraband and interband transitions, the absorption 
by magnons is two to three orders of magnitude weaker, yielding an optical conductivity of 
the order of 1 S/cm. Therefore, spin- wave excitations can be detected only in the transmit- 
tance spectra of thin single crystals. Up to now, investigations of this type have only been 
performed on single layer cuprates [|13|. Broad structures between 0.4 eV and 1.2 eV were 



observed which were interpreted as an exciton and magnon- sidebands [|13 . 

The crystals with typical dimensions of 1*1*0.1 mm^ had been annealed in the UHV at 
700 K for two days to exclude doping by excess oxygen. The samples are very close to the 
pure limit YBa2Cu306.o showing values of the conductivity function lower than 0.1 S/cm 
which is about three (five) orders of magnitude smaller than in YBa2Cu306.i (YBa2Cu307) 
in the same spectral range. The measurements were performed using a Fourier transform 
spectrometer Bruker IFS 113v in the spectral range between 85 meV and 1.5 eV. The 
samples were mounted on a diaphragm in a helium-fiow cryostat. Reference spectra at each 
temperature were obtained using a second, identical diaphragm and a turning mechanism. 
Hence an absolute photometric accuracy of the transmission data of about 1% was achieved. 

The conductivity function a{hjjj) = 2uj6on{fijuj)k{huj) can be calculated, if the sample 
thickness and both transmission and reflection spectra are known. Here, n and k denote the 
real and imaginary part of the refractive index. The sample thickness could be determined 
precisely from the spectral position of the interference maxima. Since reflection and trans- 
mission measurements use slightly different incident angles, and as furthermore the incident 
light beams are not completely parallel, there are still small interference structures in our 
plot of a^hw). 

In the upper panels of Fig. |I] reflectance and transmittance spectra, obtained at T = lOK 
on a single- crystalline platelet with a thickness of c? = 125/im, are displayed. The resulting 
conductivity function is shown in the lower panel of the same flgure. Between 0.1 eV and 
0.4 eV the measurements are dominated by interference effects, indicating regions of low 
absorption (as can be seen in a). Whereas the interference structure was precisely resolved 
by the spectrometer it cannot be resolved in the flgure. 

The reflectance spectrum by itself is not sufficient to determine the excitations present 
in this spectral range, only the knowledge of both reflectance and transmittance provides 
full information. To discuss the different excitations a plot of a{huj) is most suitable. There, 
the exponentially decreasing high energy tail of the highest fundamental phonon mode is 
observed up to 0.15 eV. Several smaller structures due to absorption by multi-phonons 
are superimposed on it. The main absorption features in the mid-infrared region are of 
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magnonic origin. The excitation of single-magnons of the optical branch is observed at 
Egp = 178.0 meV. The two peaks marked i?2a and both belong to the two-magnon 
spectrum, as will be discussed below. The broad high energy tail of the spectrum is caused 
by higher multi-magnons. Finally, the steep increase of conductivity above 1.3 eV is due to 
the onset of intrinsic absorption, the excitation of carriers across the charge transfer gap. 

In order to interpret this magnon spectrum, we use linear spin wave theory (LSW) to 
gain the excitation spectrum of localized spins on a bilayered square lattice. A Heisenberg 
Hamiltonian accounts for these low energy excitations for zero doping: 

H = J ^ ^ Sa,i^a,j + Jl2 ^ Si_iS2,i (l) 
a=l,2 <i,j> i 



where i and j label nearest neighbor sites in a two-dimensional square lattice and a G {1, 2} 
labels the two different planes in a single bilayer. Each bond is counted once. Generally 
it is found that LSW supplies quantitatively satisfying results for the Neel ground state 
at low temperatures |]T^,|T3|,|T|| . The (classical) Neel ground state Si ^ = ±(— 1)' (1, 0, 0) is 
stabihzed by a finite bilayer coupling J12 [1^ {S = 1/2). Spin-orbit effects are relatively 
small 1 16] and were neglected in Eq. (|l]). However, the finite spin-orbit coupling is needed 
to couple the external electric field to a single magnon, and further to make two-magnon 



absorption possible for the considered crystal symmetry ||T7 



Due to the finite bilayer coupling, the classical one-magnon spectrum splits into acoustic 
and optical branches [jl8l : 



huj, 



op/aci 



[k) = SJJz^ -rl + 2 (Ju/J) ■ {z ± Tfc) 



(2) 



where z=4: is the coordination number in a plane, and = 2 {cos{kxa) + cos{kya)). k is 
within the magnetic Brillouin zone, and a is the lattice constant. Absorption experiments 
probe k = with energy gap 

E,p = huopik = 0) = 2{2S)^/J^ (3) 



for the optical branch. The acoustic mode splits, due to the small spin-orbit coupling 
T6| . The gapped out-of-plane mode is indeed observed in neutron-scattering experiments 



at £'ac(fc = 0) ~4.5meV 0. The splitting of the optical branch at fc = may be estimated to 
be Aii^op ~-£'ac/2-£'op ~ l/18meV, a scale too small in comparison to the width of the optical 
magnon peak of about ImeV to be resolved in our experiment. 

The two-magnon absorption is calculated with a coupling Hamiltonian of the form 

a,h <i,j> 



where T^a,i;b,j points in the direction of the vector joining the pair (a, i; b,j) and E is the 
electric field vector |T3]. This is the only coupling allowed by crystal symmetry for a nearest 
neighbor two-magnon generation pO|. Dtz is found from a perturbation series in the long- 
wavelength electron-photon and spin-orbit interaction [O. Since we restrict ourselves to 
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fc||c, the two-magnon coupling is proportional to vrf ^ . 2 j ■ {S^ ^5'^^ — 81^82 J which creates 
a singlet pair of magnons on adjacent planes. 

Absorption is not just determined by the convoluted density of states (DOS) of two 
magnons but also by the (local) interaction between the two magnons |jl9|. Therefore, 
absorption is not just characterized by a step-like increase at the optical two-magnon edge, 
2Eop, and a diverging DOS at the upper band edge (as for J12 = 0). Rather, the local 
interaction will reduce the frustration produced by the two spinflips and allow two optical 
magnons to form a nearly bound state below 2Eop. Due to the admixture of acoustic magnons 
the bound state shows up as a resonance (see Fig. |^). A second broader resonance close to 
the band edge replaces the diverging two-magnon DOS for J12 > 0, similar to the well-known 
two-magnon peak in Big Raman scattering. Quantum fluctuations certainly broaden the 
two-magnon peaks. 

The exact positions of both peaks depend upon the ratio Ju/J, as is displayed in the 
inset of Fig. ^ Comparison with experiment yields a value of 0.55 ± 0.05 for this ratio 
(0.53 for £'2a and 0.58 for E2h)- Hence both resonances confirm the interpretation of the 
peak at E^p as a single magnon process. The calculated minimum at 2-i?op = 356.0 meV is 
also present in our measurements (see inset of Fig. ||). With i?op = 178meV and Eq. (0) we 
obtain J ^ 120 meV and Ji2~ 66meV. For comparison, the calculated absorption spectrum 
for J12/ J = 1 is given by the dashed line in Fig. ^ which demonstrates that the shape has 
changed qualitatively, quite distinct from the measured spectrum. Actually, the resonance 



at E2h is lost for J12/ J ^ 0.6 |^ and is replaced by a broad hump. 

The amazingly high value of J12 implies that the bilayer coupling may not be neglected in 
the interpretation of any experiment with YBa2Cu306. The slope of the magnon dispersion 



as measured in neutron scattering is 25*7^ Ji2/t/ — ( J-I-J12/8), but J12/8 is still a 
small correction within experimental uncertainty. Raman-scattering results should depend 
on J12 more visibly: The calculated zone boundary two-magnon peak in Big geometry is 
located at approximately Eb^^ — 2.7 J + J12 in LSW (Fig. ^). That is, the excitation is a 
nearly local two-magnon state in a single plane with energy {2Sz—l)J + 2S J12 (cf. |1^), and 
delocalization reduces the frustration in the plane. This calculated peak position differs from 
that of recent experiments |]TT| by about +15%. But, the standard Fleury-Loudon scattering 



Hamiltonian i^FL, which describes the interaction of spin pairs with light through a spin- 



exchange process, has to be drastically corrected for the high incident energies p3[. These 
are always above the charge transfer gap. Though, we only expect a minor shift of the 
two-magnon peak position. 

We propose that a major shift may result from coupling to phonons: It was shown 
that the unexpected large width AEb^^ of the Raman peak may be explained by such a 
coupling which is effective only for zone boundary magnons p^ , |25[] . This mechanism may be 
responsible for a shift of the Big Raman peak of the order of [AEb-^^J Eb^^Y but will effect 
the absorption peaks of Fig. ^ only little. 

Moreover, in Aig scattering geometry, which emphasizes the zone center, a broad reso- 
nance is observed at EA^g ~350meV |26]. We expect this zone center peak not to be shifted 



by phonons. Indeed, LSW reproduces such a broad feature with a maximum close to the 
measured Ea^^ (inset of Fig. H). To ensure that Ai^-scattering vanishes for J12 = 0, fiuctu- 
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ations of the longitudinal component SfSj have to be added to the response functions |T9 
(inset of Fig. ^, dashed line). Although a consistent calculation including these fluctuations 
is beyond LSW the trend of such a correction is clearly seen: the asymmetry of the Aig 
resonance is more pronounced, and the intensity is reduced. However, even if the shape is 
reproduced correctly, other processes, such as next nearest neighbor contributions to ifpL, 
may be as important. Besides, also single layered cuprates exhibit a Aig resonance which 
does not result from -ffpL- 

In conclusion, we presented the low conductivity mid-infrared spectrum of YBa2Cu306.o, 
which is dominated by absorption due to magnons. By comparison with linear spin wave 
theory we could determine the ratio Ju/J = 0.55 ± 0.05 and hence J ~ 120 meV and 
Ji2 ~ 66 meV. 
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FIG. 1. Upper panels: In-plane reflectance and transmittancc spectra (A;||c) of a YBa2Cu306.o 
single crystal with d = 125^m at T= 10 K. The arrows mark absorption processes described in the 
text. Lower panel: The resulting optical conductivity function. 

FIG. 2. Calculated two-magnon absorption for Ji2/J = 0.55 (solid line) and 1.0 (dashed line) 
with fc||c, r = 0. Inset: Peak positions of the lower (£?2a) and upper (£^2b) two-magnon resonance 
in units of Eop. Arrows indicate the experimental values. 

FIG. 3. Calculated Raman scattering in Big geometry for Ji2/J = 0.55, r = 0. The Raman shift 
hAuj is displayed in units of E'op to compare the peak position to the position of the absorption 
resonances in Fig. 2. Inset: Aig geometry in LSW (solid line), with corrections from longitudinal 
spin fluctuations (dashed line, scaled by a factor 3.5), both for Ji2/J = 0.55. 
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